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The  purpose  of  this  work  is  to  develop  a theoretical  analysis  of  unsteady  proces-| 
ses  in  solid  propellant  combustion,  particularly  combustion  stability  and 
extinction  by  rapid  depressurization.  It  is  assumed  that  the  solid  decomposes  at] 
its  surface  by  a pyrolysis  law,  and  the  gaseous  decomposition  products  react 
exothermically  following  an  Arrenius  law.  The  nondimensional  activation  energy 
of  the  gas  phase  reaction  is  considered  to  be  large  so  that  the  reaction  takes 
place  in  a thin  zone  where  the  temperature  is  close  to  the  flame  temperature. 

The  gas  phase  is  assumed  to  be  quasisteady  and  under  this  condition  the  equations! 
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n the  gas  phase  are  solved  in  the  limit  ot  high  activation  energies.  The  inalysis 
of  the  energy  equation  in  the  condens/ed  phase  chows -that  the  characteristic 
response  time  of  the  solid  to  gas  phase  perturbations  is  large  compared  to  the 
characteristic  residence  time  in  the  heat-up  zone  of  the  solid,  their  ratio  being 
of  the  order  the  nondimensional  activation  energy  in  the  gas  phase  reaction.  A 
stability  criterion  has  been  obtained,  resulting  that  for  large  values  of  the 
nondimensional  activation  energy  of  the  pyrolysis  law  there  is  always  unstable 
behaviour.  For  the  stable  case  the  dynamic  extinction  problem  is  being  present), 
investigated . 
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ABSTRACT 


The  purpose  of  this  work  is  to  develop  a theoretical  anal- 
ysis of  unsteady  processes  in  solid  propellant  combustion,  par- 
ticularly combustion  stability  and  extinction  by  rapid  depres- 
surization. The  interest  in  understanding  combustion  stability 
stems  from  the  fact  that  it  may  lead  to  inefficient  operation 
in  rocket  motors.  Extinction  by  depressurization  is  of  interest 
in  order  to  design  rocket  engines  for  space  application  with 
stop-restart  capability. 

It  is  assumed  that  the  solid  decomposes  at  its  surface  by 
a pyrolysis  law,  and  the  gaseous  decomposition  products  react 
exothermically  following  an  Arrhenius  lav?.  The  nondimensicnal 
activation  energy  of  the  gas  phase  reaction  is  considered  to 
be  large  so  that  the  reaction  takes  place  in  a thin  zone  where 
the  temperature  is  close  to  the  flame  temperature. 

The  gas  phase  is  assumed  to  be  quasisteady  and  under  this 
condition  the  equations  in  the  gas  phase  are  solved  in  the  limit 
of  high  activation  energies  yielding  conditions  for  the  burning 
rate  and  heat  flux  to  the  solid,  The  analysis  of  the  energy  equa- 
tion in  the  condensed  phase  shows  that  the  characteristic  res- 
ponse time  of  the  solid  to  gas  phase  perturbations  is  large  com- 
pared to  the  characteristic  residence  time  in  the  heat-up  zone 
of  the  solid,  their  ratio  being  of  the  order  of  the  nondimensio- 
nal  activation  energy  of  the  pyrolysis  law  there  is  always  uns- 
table behaviour.  For  the  stable  case  the  dynamic  extinction  pro- 
blem is  being  presently  investigated.  Preliminary  results  seem 
to  give  an  extinction  criteria  for  step-like  pressure  changes. 


Future  work  will  include  heat  losses  and  will  calculate 
higher  order  terms  in  the  high  activation  energy  limit.  Fossi- 
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NOMENCLATURE 


B Preexponential  factor 


c Specific  heat 


D Diffusion  coeficient 


E Activation  energy 


E'  Dimensionless  activation  energy 


1 Dimensionless  heat  flux  at  the  surface,  see  Eq,(14) 


L Heat  of  vaporization 


L Lewis  number 
e 


m Burning  rate 


n Exponent  of  fuel  mass  fraction  in  chemical  reaction  rate 


p Pressure 


Q Heat  released  per  unit  mass  of  fuel 


R Universal  gas  constant 


T Temperature 


t 


Time 


Y Fuel  mass  fraction 

y Dimensionless  space  coordinate 

a Thermal  difussivity 

y Dimensionless  steady  state  surface  temperature,  see  Eq,C38) 
e Dimensionless  inverse  activation  energy 
0 Nondimensional  temperature 

Thermal  conductivity 

V Dimensionless  burning  rate 

£ Nondimensional  space  coordinate 

7T  Nondimensional  pressure  history  j 

p Density 

Dimensionless  time,  see  Eq,(19) 

/ 

x Dimensionless  time,  see  Eq.(30) 

J 


See  Eq,C38) 


SUBSCRIPTS 

Flame 

Gas 

Surface 

Steady  state  conditions 


Infinity  in  the  solid 


1 . Introduction 
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An  understanding  of  the  unsteady  burning  of  solid  propellants 
provides  insight  into  such  important  problems  as  the  pressurization, 
depressurization  and  stability  of  a solid-rocket  motor. 

Research  in  the  stability  of  combustion  in  a rocket  motor, 

has  been  active  over  the  past  twenty  years,  since  irregular  pulses 

0 

in  chamber  pressure  were  observed  to  develop  instead  of  the  expected 
smooth  pressure-time  history.  These  irregular  pulses  are  generally 
accompanied  by  more  regular,  small-amplitude,  pressure  oscillations, 
with  frecuencies  of  the  order  of  the  natural  vibrational  frecuencies 
of  sound  waves  in  the  chamber.  Combustion  instability  leads  to 
inefficient  operation  of  rocket  motors  and  even  to  mechanical  failure 
of  the  propellant. 

A considerable  theoretical  effort  has  been  devoted  towards 
the  understanding  of  this  phenomenon.  However,  there  are  so  many 
different  effects  that  may  influence  the  stability  of  a burning 
solid  as  to  prevent  a consensus  on  the  theoretical  description  of 
the  phenomenon. 

An  excellent  review  of  the  analyses  of  the  small  amplitude 
. . 7 

pressure  oscillations  has  been  developed  by  Culick  . These  analyses 
calculate  the  admittance  function  of  the  burning  surface  which 
gives  the  burning  rate  response  to  a pressure  disturbance,  indicating 
therefore  wether  pressure  oscillations  are  amplified  or  attenuated. 
Most  of  the  analyses  assume  the  gas  phase  to  be  quasisteady,  in  the 
sense  that  the  gas  phase  adjusts  very  quickly  to  changes  in  condi- 
tions when  compared  to  the  response  of  the  solid,  and  differ  mainly 
in  the  assumptions  used  to  calculate  the  heat  transfer  from  the  gas 
to  the  solid  phase.  However  as  Culick  points  out,  the  majority  of 
the  results  lead  to  the  same  two-parameter  form  of  the  admittance 
function,  with  different  definitions  of  the  two  parameters.  A 
stability  boundary  in  the  space  with  coordinates  the  two  parameters 


It 


entering  the  admittance  function,  is  defined  as  the  curve  where 
the  admittance  function  becomes  infinite,  so  that  a small  pressure 
change  causes  a large  fluctuation  in  burning  rate.  Points  above 
this  stability  boundary  produce  unstable  solutions  and  points  below 
this  boundary  provide  stable  steady  solutions. 

Interest  in  the  response  of  a burning  solid  to  an  externally 
applied  pressure  variation  stems  from  the  possibility  of  extinction 

Q 

by  a rapid  pressure  decay  . This  dynamic  extinction  process  is 
useful  in  order  to  design  solid  propellant  rockets  with  stop-restart 
capabilities,  of  interest  in  connection  with  space  applications. 

Nearly  all  the  theoretical  models  developed  to  explain  dynamic 
9-13 

extinction  invoke  some  kind  of  quasisteady  approximation  for 

the  gas-phase.  However,  as  pointed  out  in  reference  (12)  some  of 
9-11 

these  analyses  have  interpreted  incorrectly  this  assumption  by 

implying  that  the  heat  feebback  from  the  flame  to  the  solid  is  a 
steady  state  function  of  the  instantaneous  pressure  only. 


The  present  paper  is  a first  step  towards  analyzing  unsteady 
processes  in  solid  propellant  burning  by  means  of  asymptotic  techni- 
ques based  on  the  assumption  that  the  nondimens ional  activation  energy 
of  the  gas-phase  reaction  is  large.  Attention  is  paid  both  to  the 
problem  of  the  stability  of  the  steady  state  and  to  the  response  of 
the  burning  propellant  to  an  externally  imposed  pressure  variation. 


We  consider  a one-dimensional  model  in  which  a condensed  material 

gasifies  by  a rate-controlled  surface  process  and  then  reacts  in  the 

gas  phase.  This  gas  phase  reaction  is  described  by  an  Arrhenius 

law,  and  we  consider  the  limit  in  which  the  nondimensional  activation 

energy  is  large.  The  quasisttady  assumption  is  used  to  describe  the 

3 

gas-phase,  so  that  we  may  use  the  results  obtained  by  Williams  and 

4 . 

Buckmaster  et  al  , when  analyzing  the  quasisteady  burning  of  a solid 

in  the  limit  of  high  activation  energy.  These  analyses  yield  the 
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burning  rate  and  the  heat  feedback  to  the  solid  as  func- 
tions of  the  pressure  and  the  flame  temperature,  and  the_ 
se  relations  are  then  used  to  analyze  the  unsteady  respon 
se  of  the  condensed  phase.  It  is  found  that  the  charac- 
teristic response  time  of  the  solid  is  large  of  the  order 
of  the  nondimensional  activation  energy,  so  that  the  temp- 
erature profiles  are  quasisteady  is  first  approximation 
and  the  transient  term  is  a perturbation  of  the  quasisteady 
solution  which  may  be  calculated  from  the  quasisteady  pro- 
files. In  this  way  a differential  equation  is  derived 
which  describes  the  evolution  with  time  of  the  burning  rate 
as  a function  of  the  nondimensional  surface  temperature, 
the  nondimensional  pyrolisis  activation  energy  and  the  pres- 
sure-time history.  The  analysis  of  this  equation  yield 
the  stability  condition  of  the  steady  state  solution,  and 
the  response  of  the  solid  to  an  imposed  pressure  variation. 
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2 . F ormulat ion 


We  consider  a one-dimensional  model  with  the  solid  occu- 
pying the  half  space  x < 0 and  the  gaseous  phase  the  region 
x > 0.  For  convenience,  the  origin  is  fixed  at  the  surface  of 
the  regressing  solid.  Figure  1 is  a schematic  representation  of 
the  process  and  shows  the  effects  accounted  for  in  this  study. 


Surface  gasification  is  assumed  to  occur  by  an  irrever- 
sible pyrolisis  process  which  is  described  by  an  Arrhenius  law. 
A one  step  over-all  exothermic  reaction  takes  place  in  the 
(premixed)  gas.  The  present  model  has  been  used  because  it  has 
been  successfull  in  describing  the  steady-state  deflagration  of 
several  propellants. 


The  equation  of  conservation  of  momentum  reduces  to  the 

statement  that  the  pressure  is  approximately  uniform  throughout 

1 

the  region  treated  but  varies  with  time  , and  we  will  also  use 
the  well  justified  assumptions  that  the  work  associated  with 
viscous  and  external  forces  is  negligible  and  we  will  use  Fick's 
law  to  calculate  the  diffusion  velocities. 


With  these  assumptions,  the  conservation  equations  of  mass, 
concentration  and  energy  in  the  condensed  phase  and  in  the  gas 
phase,  become  respectively. 
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where  all  symbols  are  defined  in  the  nomenclature 


These  equations,  together  with  the  equation  of  state  should 

determine  p , u,  Y,  T,  T . 

g g 

The  boundary  conditions  are 
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The  pyrolisis  process  is  assumed  to  follow  an  Arrhenius  law 


m = B'  exp ( -E/RT  ) (9) 

s 

where  m is  the  mass  flux  relative  to  the  burning  surface,  and  E 
the  activation  energy  of  the  pyrolisis  process. 

3 . Gas-phase  analysis 

In  most  studies  of  unsteady  solid  propellant  burning,  the 
gas-phase  is  assumed  to  be  quasisteady  in  the  sense  that  the  res- 
ponse time  in  the  gas  is  short  compared  to  the  response  time  of 
the  condensed  phase.  The  ratio  of  these  characteristic  times  is 

of  the  order  of  the  ratio  of  the  thermal  responsivities  of  solid 

2 

and  gas,  which  is  usually  small  . The  quasisteady  assumption  will 
not  be  valid  when  analyzing  the  response  of  a burning  solid  to 


I mm 


r 


10.T 


very  high  frecuency  pressure  oscillations' 


When  the  gas-phase  is  considered  to  be  quasisteady,  the  ana- 
lysis is  greatly  simplified  since  the  conservation  equations  are 
uncoupled.  All  the  time  derivate  terms  may  be  neglected  so  that 
equation  (1)  reduces  to  p u = m . In  addition  the  dp/dt  term  in 

o 

equation  (3)  may  be  neglected.  Under  these  conditions  the  gas-phase 
equations  may  be  solved  in  the  limit  of  high  activation  energy  of  the 
gas-phase  reaction.  This  quasisteady  solution  has  been  derived  by 
in  the  case  of  Lewis  number  unity  and  by  Buckmaster  et 
, for  arbitrary  Lewis  number,  Le . 


Williams 
4 


al 


Rather  than  repeat  those  analyses  we  will  only  state  the  results 
needed  in  the  following  sections.  The  reader  is  referred  to  reference 
3 and  4 for  details.  The  analysis  in  reference  (4)  considers  Le  cons- 
tant but  arbitrary.  In  the  limit  of  high  activation  energy,  it  is 
found  that  the  reaction  term  is  only  significant  in  a thin  flame 


where  T is  close  to  T^.  In  terms  cf  the  nondimen s ional  variables 


6 = c T/Q 
g g 


C = / (m/a )dx 


(10) 


the  solution  outside  of  the  flame  sheet  becomes  simply 
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Y 
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0-1  + 1 exp ( £ ) 


1 - exp{Lg(C  ~ £f)J 


for  £ < £. 


(11) 


= e 


for  £ > £, 


= 0 


(12) 


where  5^,  the  location  of  the  flame  sheet,  is  given  t>  leading 


order  by 


£f  = In  (1/1) 


(13) 
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The  parameter 


is  the  dimensionless  heat  conducted  out  of  the  gas  at  the  interface 

The  effect  of  the  solid  and  the  pyrolisis  are  concealed  in  1.  An 

overall  energy  balance  in  the  gas  phase,  provides  a relationship 

between  0 , 0C  and  1. 
s f 


1 = 1 + 0 - 0. 

s 1 


Since  the  flame  sheet  must  lie  in  the  gas,  equation  (13) 
provides  limitations  in  the  possible  values  of  1 


0 < 1 = 1 + 0 


0f  < 1 


When  1 ■>  0 the  flame  sheet  moves  to  infinity  and  no  heat  reaches 
the  condensed  phase  from  the  flame.  When  1 -»•  1 the  flame  sheet 
approaches  the  surface  and  all  the  heat  generated  at  the  flame  goes 
to  the  solid. 

Solutions  for  1 ->■  0 and  1 •>  1 are  presented  in  reference  (4) 

The  reaction  zone  is  located  in  the  vicinity  of  £f,  where 

( £ - £^)  is  of  order  0^/E'  . The  parameter  E1  is  the  nondimensional 

activation  energy  E'  = E c /RQ.  In  the  reaction  zone  the  tempera- 

ture  differs  from  the  flame  temperature  by  a small  quantity  of  order 
2 

Tj.R/E  . To  leading  order,  only  the  reactive  and  diffusive  terms 
are  important  in  this  zone.  Solution  of  the  energy  equation  with 
the  appropriate  matching  conditions  to  the  frozen  solution  outside 
of  the  flame  sheet,  provides  the  burning  rate  eigenvalue. 


In  this  way  an  expression  for  the  burning  rate  is  obtained 


12.- 


involving  the  pressure  and  the  flame  temperature,  namely 


\ no2n  + 2 /r,  ,n  •'■1  , E'.  n/2 

i=  2r(n  + l)aB6f  /E ' exp(-  )p 


(17) 


For  large  activation  energy  the  effect  of  the  exponential 
term  is  dominant  so  that  the  square  root  term  may  be  taken  as 
constant  when  analyzing  small  changes  in  the  flame  temperature  . 
Equation  (17)  coincides  with  the  Denison-Baum  formula'*. 

4 . Condensed  Phase  Analysis 

The  characteristic  time  in  the  condensed  phase  is  short 

2 

compared  with  the  characteristic  time  in  the  gas  phase  . Therefore 
the  condensed  phase  should  not  be  considered  quasisteady.  However, 
a brief  description  of  the  quasisteady  solution  will  be  presented 
before  considering  the  unsteady  analysis. 

In  terms  of  the  nondimensional  variables  foi?  the  solid 


T -T 
so  « 


2 

m c 

t = °—  t 

1 Ap 


Eq.  (4)  may  be  written 


36  36  _ 326 

3t1  P Dy  3y 7 


The  subscribt  o refers  to  the  initial,  steady  condition, 
boundary  condition  (6)  may  be  written, 


To  solve  the  quasisteady  problem,  we  may  neglect  the  time 
derivative  term  in  Eq.  (20).  The  solution  to  Eq.  (20)  is  simply 

6 = 6 exp  (uy)  (22) 

S 

so  that  from  Eq.  (21). 

cgTfo  3 cT»  + Q - L (23) 

This  equation  can  also  be  obtained  by  an  overall  energy 
balance  and  shows  that  the  flame  temperature  is  constant  in  a 
quasisteady  process.  For  a given  Tm,  Eq.  (23)  gives  the  value  of 
T^q  which  can  be  used  in  the  taison-Baum  formula,  Eq.  (17),  to 
calculate  m.  The  pyrolisis  law,  Eq.  (9), yields  T , so  that  Eqs . 
(11),  (12)  and  (22)  describe  the  complete  temperature  profile. 

Figure  (2)  shows  the  burning  rate  m under  steady  conditions; 

as  a function  of  pressure  for  a fixed  value  of  T and  therefore 

A 00 

of  TfQ.  This  curve  is  calculated  by  using  Eqs.  (17)  and  (23). 

The  surface  temperature  T is  a parameter  along  the  curves 

s 

of  figure  (2),  since  it  is  related  to  the  burning  rate  mQ  through 
the  pyrolisis  law.  However,  Tg  is  limited  by  the  inequalities  (16) 
so  that  only  a portion  of  the  curve  applies.  At  the  flame 

sheet  has  moved  off  to  infinity.  At  Pg  the  flame  sheet  has  reached 
the  surface. 
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Ws  will  now  analyze  the  evolution  with  time  of  an  initially 
steady  temperature  profile  during  an  unsteady  process.  After  some 
manipulations,  Eqs . (20)  and  (21)  provide  the  following  relations- 
hips 


U 0 = / -P-  dy 

s J -«  3 t ^ J 


T.p  -T  _ 
fo  f 

T -T 

SO  0° 


where  use  has  been  made  of  Eq.  (23).  The  left  hand  side  of  the 
preceding  equation  represents  the  difference  between  the  heat 
flux  existing  during  an  unsteady  process  and  the  one  existing  if 
the  process  was  quasisteady.  We  will  consider  unsteady  processes 
that  result  in  changes  of  order  unity  in  the  burning  rate  with 
respect  to  the  one  existing  under  quasisteady  conditions.  Equa- 
tion (17)  shows  that  in  the  limit  of  high  activation  energy  of 

the  gas  phase  reaction,  small  changes  in  the  flame  temperature  of 
2 

order  RT^/E^,  produce  variations  of  order  unity  in  the  burning 
rate.  Therefore,  the  right  hand  side  of  Eq.(24)  is  small,  so 
that  the  heat  conducted  to  the  solid  during  an  unsteady  process 
differs  by  a small  quantity  from  the  heat  conducted  if  the  process 
was  quasisteady. 

The  second  equality  of  Eq.(24)  indicates  that  the  variation 

with  time  of  the  heat  content  in  the  heat  up  zone  of  the  solid 

2 

is  small  of  order  RT_/E  . When  this  heat  content  decreases,  the 

f g 

solid  appears  from  the  gas  phase  as  a heat  source,  and  therefore 
the  flame  temperature  increases. 


' 

i 


It  is  necessary  to  wait  times  of  order  E /RTf  to  produce 

g * 

changes  of  order  unity  in  the  heat  content  of  the  heat  up  zone  in 
the  solid.  The  characteristic  response  time  of  the  solid  is 
therefore  long  compared  with  the  characteristic  residence 


15.- 


time  in  the  heat  up  zone. 


The  Denison-Baum  formula,  Eq.  (17),  yields 


T P 

_ , f.n+1  , .n/2 

U = ( — ) ( — ) exp 


fo 


2RT  fo  2RTf 


(26) 


which  in  the  limit  of  high  values  of  E„/2RTl  becomes 

g fo 


where 


t . .WW 

n GXP  { 2RT£q 


w , p \0/2 
^ - ( p — ) 

O 


Let's  define  a small  parameter  e,  as 

c 2RT2 

e = -&■ 


fo 


c E ( T -T  v 

g SO  °°  ) 

and  introduce  as  nondimens ional  time  variable 


t = ET, 


(27) 


(28) 


(29) 


(30) 


Equation  (27)  may  be  used  to  express  Eq . (20)  and  the 
boundary  condition  (24)  in  terms  of  the  new  time  variable  t,  as 


30  30  30 

e 37  + P 


3y 


(31) 


30 

3y 


= u[9s  ' e ln  ^ 


(32) 


Introducing  the  expansions 


-16.  T 


0 = 0 +60.,  + .. 

O 1 


4 = uo  + e + 


(33) 

(34) 


in  Eqs . (31)  and  (32),  expanding  and  collecting  like  powers  of  e 
the  following  equations  defining  0q  and  0.^  are  obtained 
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The  solution  to  Eq.  (35)  is  the  quasisteady  solution,  Eq . (1 
The  value  of  the  surface  temperature  may  be  written  in  terms  of  tl 
burning  rate  u,  through  the  pyrolisis  law,  Eq . (9),  resulting 


lny  + 1 


l-oilnp 


where 


Integrating  Eq.  (36)  from  - « to  the  surface,  the  following 
expression  is  derived 


dy  = 


-Veis‘Vos+Vo>  + W^ls-,,oln  T" 

/ \ \ / 


30 
o 

3 T 


(39) 
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The  quasisteady  solution,  Eq.  (22),  together  with  Eq.(37) 
may  be  ised  to  calculate  the  integral  in  the  left  hand  side  of  the 
preceding  equation,  thus  obtaining 

dp  ( Y-l ) ( 1-ulnw  )2  p3  u 

o o o , , o«  . 

= - £ 2 ln( — > (40) 

dr  Yto-Y  + l + m( y-2  )lnp  +to  In  (p  ) ir 

o o 

which  describes  the  evolution  with  time  of  the  burning  rate  during 
an  unsteady  process,  in  terms  of  two' parameters  y and  w.  For  a 
quasisteady  evolution  dpQ/dT  vanishes  and  therefore  pq  = it. 

The  value  of  the  burning  rate,  pq,  obtained  from  Eq.  (40), 
may  be  used  in  the  pyrolisis  law  to  deduce  the  surface  temperature 
history,  and  through  Eq.  (22)  the  complete  temperature  profile. 

5 . Stability  Analysis 

In  this  section  an  analysis  is  presented  of  the  stability  of 
the  steady  deflagration  of  a solid  which  undergoes  an  Arrhenius 
type  gas  phase  reaction  with  large  activation  energy. 

The  pressure  is  considered  to  be  constant,  so  that  n=l.  The 
burning  rate  of  the  steady  solution  is  mQ , so  that  initially  U0  = l» 
Let's  assume  that  at  t = 0 a perturbation  changes  the  value  of  the 
burning  rate  so  that  = 1+u'.  Eq.  (40)  shows  that 


du'  = _ _p 
dT  yw 


(41) 


Y-l 

Therefore  the  stability 
on  the  value  of  the  parameter 
than  one,  the  solution  is 
is  smaller  than  one.  In  the 


- 1 

of  the  steady  state  solution  depends 
YQ)  . If  this  oarameter  is  greater 

v — j[ 

stable,  being  unstable  when  it 
stable  case  the  perturbations  decrease 


J- 


I 

exponentially  to  zero  with  time,  while'  in  the  unstable  case  the 
perturbations  grow  exponentially.  The  parameter 


i 


i ' 


i 


Y-l 


yw 


E ( T -T 

so  00 


) 


RT 


2 

so 


A 


(42) 
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is  identical  with  the  parameter  A used  by  Denison  and  Baum  when 

analyzing  the  unstable  burning  of  solid  propellants.  These  authors 

develop  a linearized  analysis  of  the'  response  of  a burning  solid 

to  a pressure  disturbance  using  a model  which  parallels  the  one 

used  in  the  present  analysis.  The  difference  however,  is  that  in 

our  analysis  we  retain  the  nonlinear  effects  that  were  linearized 

\Q.f  their  paper 

in  reference  (5).  In  figure  (1),  Denison-Baum  show  the  stability 

boundary  which  separates  regions  of  stable  from  regions  of  unstable 

burning.  Two  parameters  A and  a,  define  the  stability  of  a given 

solution.  The  parameter  a is  essentially  identical  to  our  parameter 
that 

e,  and  ifT"  figure  (1)  we  observe  that  for  a ->  0 the  solutions  are 
stable  for  A < 1 and  unstable  for  A > 1.  Therefore  our  stability 

I 

criterion  coincides  with  the  one  of  Denison  and  Baum. 


The  instability  of  the  one  dimensional  deflagration  may  be 
clearly  interpreted,  when  the  activation  energy  of  the  pyrolisis 
is  large  ( w small).  Let's  assume  that  at  a certain  time  a pertur- 
bation causes  the  burning  rate  to  increase  (decrease)  with  respect 
to  its  steady  value.  Since  the  activation  energy  of  the  pyrolisis 
is  tirge , the  surface  temperature  will  remain  nearly  constant. 

However,  the  width  of  the  heat  up  zone  in  the  solid  will  decrease 
(increase),  so  that  the  total  heat  content  of  the  solid  will  decrease 
(increase).  The  variation  of  the  thermal  energy  of  the  solid  pro- 
duces an  increase  (decrease)  in  the  flame  temperature,  as  may  be 
seen  from  Eq . (24),  which  will  further  increase  (decrease)  the 
burning  rate.  This  self accelerating  behaviour  results  in  instability 
of  the  one-dimensional  deflagration.  When  the  activation  energy 
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of  the  pyrolisis  is  not  large,  there  are  two  effects  which  control 
the  response  of  the  burning  solid  to  a perturbation  in  the  burning 
rate.  When  the  burning  rate  increases  the  width  of  the  heat  up 
zone  decreases  and  the  surface  temperature  increases.  The  decrease 
in  width  tends  to  decrease  the  heat  content  of  the  solid  thus  pro- 
ducing instability,  and  the  increase  in  surface  temperature  tends 
to  increase  the  heat  content  of  the  solid  thus  producing  stability. 
The  stability  of  the  solution  depends  on  which  of  these  two  effects 
dominates . 

6 . Response  of  a burning  solid  to  pressure  variations. 

In  this  Section  an  analysis  is  presented  of  the  response  of 

a burning  solid  to  an  imposed  pressure  variation.  The  evolution 

with  time  of  the  burning  rate  is  given  by  the  solution  of  Eq  . (40) 

within  being  a known  function  of  time  defined  by  Eq . (28). 

the’  nondimensional  pressure 

It  was  shown  in  section  (4)  that  the  characteristic  response 
time  of  the  solid  is 


t 

c 


Xp 


2 

m c 
o g 


E (T  -T  ) 
g so  °° 


2RT 


2 

f o 


(43) 


Therefore,  if  the  characteristic  time  of  pressure  variation  is 
small  compared  with  it,  the  pressure  variation  function  t is  a 
step  function.  If  the  time  of  pressure  variation  is  long  compared 
to  the  characteristic  time  shown  in  Eq.  (43)  the  response  of  the 
solid  may  be  considered  as  a sequence  of  stationary  states,  and  at 
any  given  time  p = m . 


meters 


The  solution  to  Eq.  (40)  depends  on  the  values  of  the  para- 
Y,  a)  and  of  the  pressure  time  history. 
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To  analyze  the  different  types  of  solutions  that  may  result 
from  integration  of  Eq.  (40),  let's  first  study  the  singularities 
of  that  equation. 

The  numerator  vanishes  at 


U = exp  ( 1 /<j)  ) 


y = it 


V = 0 


and  the  denominator  at 


In  y = - 


{ Y ~ 2 + /yCy-'+w)  } 


In  the  present  analysis  we  will  only  consider  depressurization 
problems,  so  that  it  will  decrease  from  1 at  time  0 to  its  final 
value  ti  . In  addition  we  will  always  assume  that  at  the  initial 
pressure  the  solution  is  stable,  so  that  yw/(y-l)  is  greater  than 
one.  The  existence  of  singularities  (where  the  denominator  and 
numerator  vanish  simultaneously)  depends  on  the  values  of  Y » w 
and  and  it  is  observed  that  in  order  to  have  real  solutions  to 

Eq.  (45)  it  is  necessary  that 

Y > 4to  ( 46  ) 

In  dimensional  variables  E /RT  >1  for  the  existence  of  roots. 
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It  is  observed  that  when  the  solution  at  the  initial  pressure 
is  stable,  both  roots  of  Eq.(45)  have  the  same  sign.  So  that 


Y < 2 

Y > 2 


2 positive  roots'^ 
2 negative  roots  J" 


and  the  denominator  vanishes  at  two  different  values  of  y,  which 
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are  either  both  greater  than  one  or  both  smaller  than  one.  In 
dimensional  variables  Tgo  > 2Tra  for  the  existence  of  two  negative 
roots  . 


Since  n is  a decreasing  function  of  time,  which  starts  at  1 
and  goes  to  rr^.,  singularities  will  only  exist  when  y > 4oj,  y > 2 
and  smaller  than  the  larger  root  of  Eq . (45).  When  rrf  is  smaller 
than  both  roots  there  are  two  singularities,  and  when  tr ^ lies  between 
both  roots  there  is  only  one  singularity . 

The  condition  for  the  steady  state  solution  at  the  final 
pressure  being  stable,  may  be  written 

Yu 

> i (48) 

(1-wln  ir^My-l  + w In  ir^) 

The  left  hand  side  of  this  equation  is  represented  in  Figure 
. When  In  nf  . 

(3)  as  a function  of  In  it ^ .approaches  the  value  1/w  (T  -*  <*>)  the 

left  hand  side  of  Eq.  (48)  approaches  °°  , and  the  solution  at  the 

final  pressure  is  stable.  Analogously  the  solution  is  stable  when 

In  t\ f approaches  the  value  (l-y)/w  (Tg  -*■  0).  Between  these  two 

values,  the  function  has  a minimum  when  lnu^  equals  (2-y)/2ui  and 

depending  on  wether  the  value  of  the  function  at  this  minimum  is 

larger  or  smaller  than  one,  the  solution  will  be  stable  for  any 

final  pressure  or  there  will  be  a range  of  final  pressures  for 

which  the  solution  will  be  unstable.  The  equation  which  gives  the 

values  of  In  tt^.  where  the  function  defined  in  Eq . (48)  crosses  to 

the  unstable  region  coincides  with  Eq.  (45). 

Let's  denote  these  two  values  (if  they  exist)  by  In  and 
In  it  ^ - We  may  conclude  that  if  y < 4w  the  solution  at  any  final 
pressure  is  stable.  Also,  if  y > 4i»)  and  y < 2,  the  solution  is 
stable  for  any  final  pressure  smaller  than  the  initial  pressure. 


being  unstable  in  a finite  range  of  pressures  larger  than  the  initial 
pressure.  Finally,  if  v > Hu  and  y > 2 the  solution  is  unstable  in 
a range  of  values  of  final  pressures  smaller  than  the  initial  pres-  I 
sure.  We  may  briefly  summarize  the  different  possible  cases  in  a ’] 
depressurization  process  during  which  the  nondimensional  pressure  li 
goes  from  the  initial  value  of  1 , to  a final  value  it  < 1.  j 


a)  y < 4w (- 


> 1 ) 


E : The  steady  state  solution  is  stable  at  the 

initial,  final  and  all  intermediate  pressures  occurring  during 
the  process.  Figure  (4)  shows  the  evolution  with  time  of  the 
burning  rate  for  the  case  in  which  a step  change  in  pressure 
takes  place  at  time  zero.  Figure  (5)  shows  the  burning  rate 
evolution  for  depressurization  processes  of  the  type 

“ Id  t 

ir=  (l-a)e  + a (4g) 

with  different  values  of  b. 


b)  y>4(o,y<2(- 


< l T < 2T  ) 

so °° • : This  case  is  analogous 


to  the  preceding  one.  The  steady  state  solution  at  the  final 
pressure  is  more  stable  then  the  steady  state  solution  at  the 
initial  pressure,  so  that  a perturbation  dies  out  more  rapidly 
at  the  final  than  at  the  initial  pressure. 

c)  Y > 4(o,  y > 2,n^>  ir^  : This  case  is  analogous  to  case  (a)  with 
the  solution  at  the  initial  pressure  being  more  stable  than  the 
solution  at  the  final  pressure. 


d)  Y > 4(0, y > 2,w1>  irf>  tt2> 


The  solution  is  unstable  at  the  final 


pressure,  and  Eq.  (40)  has  a singularity  at  p = tt  , and 


such  that  t(t^)  = The  character  of  the  singularity  depends 

on  the  value  of  the  slope  of  the  pressure  time  function  curve 


a 


r 
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dir 


dT 


( T = T„  ) 


Y-1  p^(l-wlnp  )2 


/y  (y-4ii)  ) 


(50) 


the  singularity  is  a node,  being  a spiral  when  this  inequality  is 
not  satisfied.  Figure  (6)  shows  the  isoclines  of  Eq.  (40)  for  the 
two  possible  cases,  together  with  the  solution  with  the  initial  con- 
dition p = 1.  When  the  singularity  is  a node,  the  solution  to  Eq. 

(40)  will  be  a smooth  function  approaching  tt,.  for  t -*■  « . However 

\£qt 

it  is  observed  that  p < the  solution  enters  a region  of  unstable 
states  and  any  perturbation  will  grow  exponentially  with  time.  The 
one-dimensional  flame  is  unstable,  so  that  our  solution  is  no  longer 
valid.  One  may  speculate  that  when  the  flame  becomes  unstable,  some 
type  of  two-dimensional  effect  or  other  effects,  will  appear  in 
order  to  stabilize  the  flame.  When  the  singularity  is  a spiral,  our 
solution  ceases  to  be  valid  when  p approaches  , since  the  variation 
of  burning  rate  with  time  becomes  very  rapid,  so  that  it  is  no  longer 
valid  to  assume  that  the  time  derivative  term  is  only  a perturbation 
of  the  quasisteady  solution.  A different  analysis  will  be  required 
to  study  the  behaviour  of  the  solution  near  p = it^.  However,  since  p 
is  decreasing  very  rapidly  to  zero,  this  condition  may  be  identified 
as  a dynamic  extinction  caused  by  a rapid  depressurization  process 


e)  y>4w,y>2 


*f  < 


The  solution  is  stable  at  the  initial 


and  at  the  final  pressure,  but  there  is  an  intermediate  range  of 
pressures( < it<  ) where  the  solution  is  unstable.  There  are 
two  singularities  of  Eq.  (40)  at  p = n^,  t>(t^)  = tt^  and  at  4 = ^2  > 

w(t  2^  = ri  2’  T^e  singularity  at  t ^ is  a node  when  Eq  . ( 50  ) is 
satisfied  and  a spiral  in  all  other  cases.  The  singularity  at 
t j is  always  a saddle.  Figure  (7)  shows  tne  isoclines  for  the 
two  possible  cases,  and  the  solution  of  Eq.  (40) 
at  the  node,  are  given  by 


The  directions 
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where 


2u)/y(y-4u))~ 

(Y-1  )p^(l-iolnp1  )2 


(52) 


w.-'th  the  minus  sign  in  Eq.  (51)  corresponding  to  the  principal  direc- 
tion. It  is  observed  that  the  principal  direction  has  a slope  larger 
in  absolute  value  than  the  slope  of  the  pressure  function  it.  The 
directions  at  the  saddle  are  given  by 


dp  _ 
dT 


■nr 


2M 


HT 

dx 


(53) 


where 

M = 2»A(  (54) 

(y-1  )p2  ( l-tolnp2  ) 

so  that  one  direction  has  positive  slope  and  the  other  one  negative 
slope.  When  the  singularity  at  t ^ is  a node,  the  solution  starting 
at  p = 1 is  the  separatrix.  This  solution  enters  the  unstable  re- 
gion through  the  node  and  exits  by  the  saddle,  reaching  in  this  way 
the  region  of  stable  states.  However  any  perturbation  occurring 
when  the  value  of  p lies  in  the  unstable  region  will  grow  with  time 
so  that  the  one-dimensional  flame  will  be  unstable  in  this  zone. 

Other  effects  not  considered  in  this  study  may  then  appear  to  stabilize 
the  flame.  As  in  the  previous  case,  if  the  singularity  at  t ^ is  a 
spiral,  our  analysis  ceases  to  be  valid  when  p approaches  n . Tbt 
burning  rate  decreases  very  rapidly  with  time,  suggesting  the  onset 
of  dynamic  extinction. 
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Figure  (8)  shows  the  stability  regions  as  a function  of  the 
parameters  y and  w,  At  the  initial  point  of  a depressurization 
process  (in  pc 0),  the  steady  state  solution  is  stable  or  not 
depending  on  wether  the  value  of  w is  at  the  right  or  left  of 
the  intersection  point  of  the  y curve  with  the  w axis  respec- 
tively. The  steady  state  solution  at  the  final  pressure  of  a 
depressurization  process,  has  the  same  value  of  y and  w but 
a different  value  of  p(ln  p<0).  The ' stability  of  this  final 
point  depends  on  wether  or  not  it  is  located  at  the  right  or 
left  of  the  corresponding  y curve.  It  is  observed  that  for 
y>2  there  is  a vertical  tangent  of  the  y curves,  and  therefore 
<’  it  is  possible  to  have  depressurization  processes  such  that 


the  steady  state  solutions  corresponding  to  the  initial  and 
final  pressure  are  stable  although  there  are  intermediate  pres- 


7.  Conclusions 


We  have  developed  an  asymptotic  analysis  of  the  unsteady  burn- 
ing of  a solid  propellant  in  the  limit  of  high  activation  energy 
The  analysis  shows  that  the  ratio  of  the  characteristic  pressure 
and  thermal  response  times  of  the  solid,  (see  Eq.  (43))  is  of  the 
order  of  the  nondimens ional  activation  energy  of  the  gas  phase 
reaction.  By  integrating  the  condensed  phase  energy  equation,  an 
equation  is  derived  whose  solution  yields  the  evolution  with  time 
of  the  burning  rate. 

2 

It  is  found  that  for  E ( T -T  )/RT  greater  than  one  the  steady 

so  °°  so 

• • • 1 

state  solution  is  dinamically  unstable.  However,  it  is  not  clear 

1 what  happens  when  the  one -dimens ional , steady  deflagration  solution 

1 4 

becomes  unstable.  T'ien  , identifies  the  onset  of  instabilities 
with  the  occurrence  of  extinction,  thus  defining  a low  pressure 
' dynamic  deflagration  limit,  which  does  not  coincide  with  the  low 

pressure  static  deflagration  limit.  It  is  more  likely  however,  that 
some  nonlinear  mechanism  or  tri-dimensional  effect  will  prevent  the 
growth  of  the  oscillations  and  the  occurrence  of  extinction,  so  that 

I 

only  the  cne-dimensional  steady  deflagration  described  in  the  present 
* analysis  will  cease  to  exist,  not  the  flame.  The  intrinsic  instabili- 

ty of  the  one-dimensional  deflagration  is  analogous  to  the  problem 

of  normal  flame  propagation  which  has  been  recently  analyzed  by 
15 

Sivashinsky  , showing  that  for  large  activation  energy  and  Lewis 
number  greater  than  one  the  solution  becomes  unstable. 

The  analysis  has  also  been  applied  to  study  the  response  of  a 
burning  solid  subjected  to  a depressurization  process.  The  character- 
istic response  time  of  the  solid  is  given  be  Eq.(43),  so  that  when 
the  time  of  pressure  variation  is  short  compared  with  it,  the  de- 
pressurization process  appears  as  a step  change  in  pressure.  We 
assume  that  the  steady  solution  at  the  initial  pressure  is  stable, 
and  derive  the  conditions  necessary  for  the  steady  solution  to  be 
stable  at  the  final  pressure  and  at  all  intermediate  pressures. 


It 
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is  shown  that  when  the  steady  solution  is  stable  at  all  pressures 

occuring  during  the  process,  the  burning  rate  always  approaches 

for  longe  times,  the  value  corresponding  to  the  steady  solution  at 

the  final  pressure.  The  burning  rate  does  not  adjust  instantaneously 

to  the  applied  pressure  so  that  during  a fast  change  in  pressure 

the  value  of  the  burning  rate  is  larger  that  the  quasisteady  value, 

approaching  it  for  long  time.  The  absence  of  dynamic  extinction 

1 3 

appears  to  coincide  with  the  suggestion  of  T'ien  , who  considers 
that  heat  losses  is  the  mechanism  responsible  for  dynamic  extinction. 

From  analysis  of  Eq.(40)  which  gives  the  time  evolution  of 
the  burning  rate,  it  is  found  that  when  the  steady  solution  is 
stable  at  the  initial  pressure  but  becomes  unstable  at  some  lower 
pressure  which  is  reached  during  the  process,  a singularity  of  Eq. 
(40)  occurs  when  the  pressure  reaches  the  value  at  which  the  steady 
solution  becomes  unstable.  The  character  of  the  singularity  depends 
on  the  slope  of  the  pressure-time  curve  at  the  singularity.  When 
the  slope  is  larger  than  a critical  value  defined  in  Eq.(50)  the 
singularity  is  a spiral.  The  solution  derived  in  the  present  analysil 
is  valid  until  a certain  time  when  the  variation  of  the  burning 
rate  with  time  becomes  infinite.  From  that  time  the  transient  term 
becomes  of  the  order  of  the  convective  and  diffusive  terms  and  our 
analysis  ceases  to  be  valid.  However,  since  the  burning  rate  is 
rapidly  decreasing  towards  zero,  this  condition  may  be  identified 
with  the  onset  of  dynamic  extinction.  When  the  slope  of  the  pressure* 
time  curve  is  smaller  than  the  critical  value  defined  in  Eq.(50), 
the  singularity  is  a node.  The  burning  rate  solution  in  absence  of 
disturbances  will  approach  for  long  times  its  steady  value  corre- 
sponding to  the  final  pressure.  However,  in  the  region  of  unstable 
states,  an‘y  disturbance  will  be  amplified  so  that  the  steady  one- 
dimensional deflagration  ceases  to  be  a valid  solution,  and  some 
other  type  of  solul^ion  will  exist. 


r 


REFERENCES 


1.  Williams,  F.A.,  "Response  of  a Burning  Solid  to  Small-Amplitude 
Pressure  Oscillations",  J.  Appl.  Phys  . , 33,  3153-3166  (1962  ). 

2.  Summerfield,  M.,  Caveny,  L.H.,  Battista,  R.A.,  Kubota,  N., 
Gostintsev,  Yu.  A.  and  Isoda,  H.,  "Theory  of  Dynamic  Extinguish 
ment  of  Solid  Propellants  with  Special  Reference  to  Nonsteady 
Heat  Feedback  Law",  AIAA  paper  70-667  (1970). 

3.  Williams,  F.A.  "Quasi-Steady  Gas-Phase  Flame  Theory  in  Unsteady 
Burning  of  a Homogeneous  Solid  Propellant",  AIAA  J.,  11, 
1328-1330  (1973). 

4.  Buckmaster,  J.D.,  Kapila,  A.K.  and  Ludford,  G.S.S.,  "Linear 
Condensate  Deflagration  for  Large  Activation  Energy",  Acta 
Astronaut ica , 3,  593-614  (1976). 

5.  Denison,  M.R.  and  Baum,  E.,  "A  Simplified  Model  of  Unstable 
Burning  in  Solid  Propellants",  ARS  Journal,  31,  1112-1122  (1961) 

6.  Green,  L.  Jr.,  Jet  Propulsion,  26,  655  (1956). 

7.  Culick,  F.E.C.,  "A  Review  of  Calculations  for  Unsteady  Burning 
of  a Solid  Propellant",  AIAA  J.,  6,  2241-2255  C1968). 

8.  Ciepluch,  C.C.,  "Effect  of  Rapid  Pressure  Decay  on  Solid  Propel- 
lant Combustion",  ARS  Journal,  31,  1584-1586  (1961). 

9.  Von  Elbe,  G.,  "Theory  of  Solid  Propellant  Ignition  and  Response 
.-■to  Pressure  Transients",  Interagency  Solid  Propulsion  Meeting 

3,  95  (1963). 


Paul,  B.E.,  Lovine,  R.L.  and  Fong,  L.Y.,  "A  Ballistic  Explana- 
tion of  the  Ignition  Pressure  Peak",  AIAA  paper  64-121  (1964). 

Horton,  M.D.,  Bruno,  P.S.  and  Graesser,  E.C.,  "Depressurization 
Induced  Extinction  of  Burning  Solid  Propellant",  AIAA  J.,  6, 
292-297  (1968). 

Merkle,  C.L.,  Turk,  S.L.  and  Summerfield,  M.,  "Extinguishment 
of  Solid  Propellants  by  Depressurization:  Effects  of . Propellant 
Parameters",  AIAA  paper  69-176  (1969). 

T'ien,  J.S.,  "A  Theoretical  Criterion  for  Dynamic  Extinction 
of  Solid  Propellants  by  Fast-Depressurization",  Comb.  Sci.  6 
Tech.,  9,  37-39  (1974). 

Baliga,  B.R.  and  T'ien,  J.S.,  "A  Theory  of  Flammability  Limits 
of  Solid  Propellants  at  Low  Pressure",  AIAA  paper  75-229  (1975) 

Sivashinsky,  G.I.,  "On  a Converging  Spherical  Flame  Front", 

Int.  J.  Heat  6 Mass  Transfer,  17,  1499-1506  (1974). 


Figure  6b  Isoclines  of  equat ion ( 40 ) when  there  is  only  a singu- 
j larity  which  is  a spiral 

Figure  7a  Isoclines  of  equation(40)  when  there  are  two  singu- 
larities and  the  upper  one  is  a node 

Figure  7b  Isoclines  of  equat ion (X0 } when  there  ape  two  singur 
larities  and  the  upper  one  is  a spiral 

Figure  8 Stability  regions  as  a function  of  the  parameters 

Y and  u> 


